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ABSTRACT 

Binarity is believed to dramatically affect the history and geometry of mass 
loss in AGB and post- AGB stars, but observational evidence of binarity is sorely 
lacking. As part of a project to look for hot binary companions to cool AGB 
stars using the GALEX archive, we have discovered a late-M star, YGem, to 
be a source of strong and variable UV emission. YGem is a prime example of 
the success of our technique of UV imaging of AGB stars in order to search for 
binary companions. Y Gem's large and variable UV flux makes it one of the most 
prominent examples of a late AGB star with a mass accreting binary companion. 
The UV emission is most likely due to emission associated with accretion activity 
and a disk around a main-sequence companion star. The physical mechanism 
generating the UV emission is extremely energetic, with an integrated luminosity 
of a fewxL at its peak. We also find weak CO J=2-l emission from Y Gem with 
a very narrow line profile (FWHM of 3.4km s _1 ). Such a narrow line is unlikely 
to arise in an outflow, and is consistent with emission from an orbiting, molecular 
reservoir of radius 300 AU. Y Gem may be the progenitor of the class of post- AGB 
stars which are binaries and possess disks but no outflows. 

Subject headings: binaries: general, stars: AGB and post-AGB, stars: mass-loss, 
stars: individual (YGem), circumstellar matter 

1. Introduction 

Binarity, very common amongst pre-main-sequence stars (e.g., Bodenheimer et al. 2000), 
is widely believed to dramatically affect the history and geometry of mass loss in late- AGB 
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and post-AGB stars (e.g., Balick & Frank 2002). The evolutionary transition from the AGB 
to the planetary nebula (PN) phase is accompanied by significant changes in the morphology 
of these objects - the roughly round circumstellar mass-loss envelopes (CSEs) of AGB stars 
evolve into pre-planetary-nebulae (PPNs) and PNs with a dazzling variety of shapes and 
intriguing aspherical symmetries (e.g. Schwarz et al. 1992, Sahai & Trauger 1998, Sahai et 
al. 2007, Sahai et al. 2011). Binarity provides a source of angular momentum, as well as 
a preferred axis to a stellar system. Critical reviews (e.g., Soker 1998) of the properties 
of bipolar PNs lead to the conclusion that binary models can explain all these properties, 
whereas single-star models have many difficulties. 

However, direct observational evidence for binarity in dying intermediate-mass stars is 
difficult to obtain. AGB stars are very luminous (~10 3 -10 4 L Q ) and surrounded by dusty 
envelopes, whereas nearby stellar companions are generally likely to be significantly less 
luminous main-sequence stars or white dwarfs. Indirect techniques such as radial-velocity 
(RV) measurements or photometric-variability (PV) measurements cannot be used for AGB 
stars since their strong periodic or irregular variability intrinsic to their pulsating atmo- 
spheres masks the corresponding variability due to a companion. PV measurements of the 
central stars of PNs (CSPNs) have resulted in detecting a total of about 40 binaries (e.g., 
Bond 2000, Miszalski et al. 2011), implying a 10-15% fraction of detectable close binaries 
among randomly selected PNs. RV measurements of CSPNs and post-AGB (pAGB) stars 
with prominent outflows (i.e., PPNs) have had very limited success (e.g., de Marco 2009, 
Hrivnak et al. 2011) because of the unknown contribution from stellar variability. Only 
for disk-prominent pAGB stars (dp AGB: Sahai et al. 2011), RV techniques have reached a 
100% success rate in detecting binarity (e.g., Van Winckel et al. 2009). 

Thus new techniques to search for the presence of binary companions in AGB stars 
are of fundamental importance. Currently, ultraviolet photometry appears to be the most 
promising technique of being able to discover a substantial number of binary AGB stars. 
An as yet untested method which could become important is to look for X-rays from AGB 
companions which have been spun- up via accretion (Soker & Kastner 2002). Most AGB 
stars with substantial mass- loss are relatively cool {T e ff < 3000K) objects (spectral types 
~M6 or later), whereas any stellar companions and/or accretion disks around them are likely 
significantly hotter (T e jj > 6000K). Observed and model spectra of AGB stars show that 
their fluxes die rapidly at wavelengths shortwards of about 2800A, hence significantly favor- 
able secondary/accretion disk-to-primary flux contrast ratios (>10) are reached in the NUV 
window of GALEX for companions of spectral type hotter than about GO (T e jj=6000K) 
and luminosity, L > 1L Q . We used UV imaging with GALEX to carry out a pilot search for 
binary companions in AGB stars (Sahai et al. 2008; hereafter Setal08). In this survey, we 
detected all of our 21 targets in the GALEX NUV filter with signal-to- noise ratios (SNR) 
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>20; 9 sources were detected in the FUV band with SNR> 8a. In this Letter, we report the 
discovery of UV emission from the AGB star, Y Gem, from archival GALEX data. 

2. Observations & Results 

Y Gem was in the field-of-view of the GALEX Medium Imaging Survey (MIS) and All- 
Sky Imaging Survey (AIS) during 3 epochs from Jan 2006 to Jan 2008 that are roughly 
spaced apart by 1 year. We discovered Y Gem's remarkable UV emission serendipitously 
while examining the GALEX archive for UV-emitting AGB stars. YGem, otherwise an 
unremarkable, semi-regular pulsating star of spectral type M8, turns out to be the brightest 
UV-emitting object amongst our sample of about one hundred M4-M8 stars with detected 
FUV emission (Sahai et al., in preparation). Following this discovery, we searched for CO 
J=2-l line emission from YGem on 2011 February 11 and 15, using the ARO 10-m SMT 
on Mt. Graham, AZ, with the 1 mm dual-polarization receiver employing ALMA Band 6 
SBS mixers. Typical system temperatures were 170-250 K. The beam size was 9 b = 32", 
and pointing accuracy is estimated to be about ±4". Observations were conducted in beam- 
switching mode with a ±3'subreflector throw. 

Y Gem was detected in both UV bands in each epoch, but its ultraviolet flux decreased 
dramatically from 2006 to 2008 (Table 1) - i.e., by a factor 12 (5.6) in the FUV(NUV) 
(Figure la). A model atmosphere corresponding to an M8III spectral type from Lejeune 
et al. (1997), fitted to the spectral-energy-distribution (SED) from blue to the far-infrared 
shows that the UV flux densities are much larger than expected from Y Gem's photosphere 
(Figure lb). Several of the flux values in 2006 & 2007 are large enough that non-linearity 
corrections, as estimated from Morrissey et al. (2007), were significant and were used to 
correct the former; the corrected values are also provided in Table 1. An examination of 
the time-tagged photon lists (averaged over 10s bins), showed no significant variability of 
the UV flux densities within each epoch, i.e., over periods of about 110 s during Epoch 1 
and 2, and 1531 s during Epoch 3. The FUV flux appears to decrease much faster than 
the NUV, with the ratio of the FUV-to-NUV wavelength-specific fluxes, Rfuv/nuv, taking 
uncorrected (corrected) values of 5.53 (4.54), 1.38 (1.47), and 2.15 (2.15) in 2006, 2007 and 
2008, respectively. Y Gem had a FUV flux of 16 mJy in 2006, much above the range of fluxes, 
~ 0.01 — 0.1 mJy, observed in the 9/21 sources in our pilot survey with FUV detections. 

The value of Rfuv/nuv can be used to estimate the equivalent blackbody temperature 
of the emitting region, T&&. Since it is difficult to assess the relative line and continuum 
contributions to the FUV and NUV fluxes, we follow Robinson et al.'s (2005) arguments for 
arbitrarily (but conservatively) assuming that half (all) of the FUV (NUV) flux is continuum 
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emission. Thus, we find that in 2006, T bb = (3 — 3.8) x 10 4 K (for the uncorrected-to-corrected 
range of Rfuv/nuv values) when the UV flux was at its peak, and T bb = 1.7 x 10 4 K, when the 
UV flux was at its lowest. The physical mechanism generating the UV emission is extremely 
energetic, with an integrated luminosity of (2.4 — 4)L at its peak in 2006. 

We detected weak CO J=2-l emission from YGem (Figure 2). The line profile appears 
centrally peaked, with a full-width at zero intensity (FWZI) of 5 km s _1 and an integrated flux 
of 0.044Kkms (1.36Jykm s _1 assuming an unresolved source). The CO line-width (FWZI) 
in YGem is quite small, compared to typical values for dusty outflows in AGB stars, >20 
km s _1 (e.g., CO catalog by Loup et al. 1993). Although we cannot rule out the possibility 
that Y Gem's CO emission comes from an anomalously slow outflow, it is more likely that 
it arises in a large, orbiting reservoir of molecular material in this object. Jura & Kahane 
(1999) argue that such narrow lines, also seen in a few other evolved stars (e.g., EUAnd), 
most likely arise from long-lived reservoirs of orbiting gas of size (100-1000) AU. The dpAGB 
star, AC Her, a known binary with a compact disk (Gielen et al. 2007), also shows a similarly 
narrow CO line (Alcolea & Bujarrabal 1991). 

We estimate a distance D=580 pc to Y Gem from its TMSS K-band magnitude m K = 
1.22, following Kahane & Jura (1994) who assume that late-M semi- regular stars have abso- 
lute magnitudes Mk = —7.6. Assuming that the gas resides in a Keplerian disk around a star 
of IMq, the speed of 1.7 km s^ 1 (assumed to be half the FWHM of the CO line) indicates an 
orbital radius of 300 AU. At this distance from the star (with luminosity of 5800L Q ), a black- 
body grain would be at a temperature of = 140 K. We adopt this temperature for the cir- 
cumstellar gas, following Jura & Kahane's (1999) similar assumption for the disk in EUAnd. 
Assuming optically-thin emission, the total gas mass is M g = 2 x 10~ 5 (4 x 10~ 4 / f C o)M Q , 
where f C o is the fractional CO-to-H 2 abundance. Given this low gas mass, the fact that 
we can discern no dust emission in Y Gem's SED - the mid- and far-infared fluxes are well 
fitted with photospheric emission from the M 8 primary (Figure lb) - is not surprising. As- 
suming a typical circumstellar gas-to-dust ratio for AGB stars of 200, the expected 60 and 
90 /im fluxes from 140 K dust are 0.5 and 0.1 Jy, significantly smaller than the photospheric 
flux at these wavelengths. 



3. Discussion 

Noting that the observed FUV fluxes in their sample were at least a factor ~30 higher 
than that can be produced by the continuum of the cool primary AGB star, Setal08 presented 
two plausible models for the FUV excesses - (1) the excess is due to the photosphere of a 
hot companion star, and (2) the excess is due to line and/or continuum emission associated 
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with accretion activity and a disk around a companion star - the accretion shock may reside 
on (a) the disk or (b) the stellar surface of the companion (as, e.g., in T Tauri stars: Calvet 
& Gullbring 1998). A 3rd model is emission from a companion related to either (a) flare 
activity on an M dwarf (dM) companion, or (b) coronal (magnetic) activity of a late-type 
(most likely ~F7-M4) main-sequence companion which has been spun-up by wind accretion 
from the primary AGB star (Soker & Kastner 2002). A 4th model is variable chromospheric- 
type line emission from the AGB star as has been inferred in the C star, TW Hor, based on 
its IUE spectra (Querci & Querci 1985). A 5th model follows from the proposal by Soker & 
Kastner (2003) of the possibility of long-duration flares in single AGB stars due to magnetic 
reconnection events leading to X-ray emission. The only AGB star with an observed X-ray 
transient is Mira A which is in a (symbiotic) binary system (Karovska et al. 2005); and 
a sensitive search for X-rays in two (apparently) single AGB stars yielded upper limits at 
levels less than (1-10)% of expectations (Kastner & Soker 2004. Hence it is not clear whether 
magnetic reconnection events in single AGB stars can produce X-ray flares; however, it would 
be worth investigating whether the latter can pwer the kind of strong, variable UV emission 
seen in Y Gem. 

We think model 3a is unlikely. Flares on M dwarfs have been seen in the UV using 
GALEX (Welsh et al. 2007), however the decay time-scales for the flare emission tends to 
be rather short - e.g., Fig. 2 of Welsh et al. (2007) shows that most M-dwarf flare sources 
decay to half their peak flux in <100s. In contrast, Y Gem's flux showed no significant 
variations during any of the 3 epochs. Furthermore, dM flares are significantly less luminous 
in the NUV and FUV than YGem: e.g., scaled to Y Gem's estimated distance of 580 pc (see 
below), the peak NUV flux of the most luminous flare object (J023955.52 — 072855.4) out 
of a total of 49 dM stars in Welsh et al's survey (who only report NUV fluxes) is roughly a 
factor 35 fainter than Y Gem's peak (corrected) NUV flux, and the giant flare in the dMe 
star G J 3685 A is a factor 80 fainter in the FUV than YGem (Robinson et al. 2005). 

Model 3b was proposed by Soker & Kastner (2002) in order to explain the X-ray emission 
observed from the central binary stars of PNs and one can imagine extrapolating it to predict 
variable chromospheric UV emission from the spun-up companion as well. However, we 
do not think models 3b and 4 are plausible models for YGem because for chromospheric 
emission, the ratio of the FUV-to-NUV fluxes (in mJy units) as measured in the GALEX 
bands, rpuv/Nuv, is much smaller than unity (about 0.1-0.05 in TWHor, and about 0.1 
in the red supergiant aOri), and results from the presence of strong emission lines in the 
NUV band. The large value of r FUV / NUV in YGem cannot be a result of dust attenuation 
because this would either not significantly change the FUV-NUV color, or would redden the 
UV spectrum. 
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The highly variable UV flux in Y Gem supports model 2 above for this source; in this 
scenario, the UV emission from YGem may indicate episodic accretion (a) onto the disk, 
from matter ejected by the AGB star (e.g., Mastrodemos & Morris 1998), or (b) onto the 
companion from the inner disk region. A likely example of the type 2a mechanis is provided 
by Mira, in which the compact companion (separation ~60AU) accretes matter in a disk 
from the primary's wind. Mira shows variable NUV and FUV line emission: UV lines 
observed with IUE during 1979-95 faded by a factor >20 by 1999-2001, and then started 
increasing back to their original levels by 2004 (Wood & Karovska 2006); the continuum was 
also variable. These authors concluded that the UV variability results from variations in 
the accretion rate (Mira A contributes to the UV flux only at wavelengths >2600A and only 
near its maximum). Mira was observed with GALEX on two different epochs separated by 
3 years, and we find that its FUV and NUV fluxes changed dramatically over this period as 
well (Table 1). 

Noting that YGem has never been classified as a symbiotic star, it is likely that its 
companion is a main-sequence star. If the emission is of type 2b, this inference is supported 
by our estimate of the emitting region's area, A emm , which is 3.3 x 10~ 3 and (2.1 — 3.3) x 10~ 3 
times the total surface area of a Sun-sized companion in 2008 and 2006, respectively (based 
on its UV flux and inferred blackbody temperature) - it exceeds the surface area of a white 
dwarf by at least a factor 5. We note that, in published (magnetospheric) accretion- shock 
models for TTauri stars (Gullbring et al. 2000, Calvet & Gullbring 1998), even the largest 
values of the energy flux of the accretion flow, T ~ 10 12 ergcm~ 2 s _1 (which determines the 
spectral shape of the UV emission), cannot produce the lowest value of Rfuv/nuv seen m 
Y Gem. Thus, if Y Gem's UV emission is due to a magnetospheric accretion- shock on the 
companion, then J 7 , and thus the mass accretion rate - typically (0.1 — few) x 10 _7 M Q in 
these models - has to be significantly higher than in T Tauri stars. 

Y Gem's relatively blue color over the wavelength region spanning the GALEX NUV- 
FUV bands is very similar to that of the carbon star V Hya, in which Setal08 found tfuv/nuv ~ 
1.1 in two separate epochs. VHya is probably the best example to-date of an evolved star 
with an active, collimated outflow, dense equatorially-flattened structures possibly related 
to a central accretion disk, and an inferred binary companion and/or accretion disk from 
UV excesses (Setal08) 

YGem appears to be one of the most prominent examples of a late AGB star with a 
(inferred) mass-accreting binary companion, but is not alone, given the results of our survey 
(Setal08). However, being the brightest, it is the best one to begin investigating in order to 
ultimately determine whether the kind of binary interaction observed in Mira is generalized 
among UV-variable AGB stars - such binary interactions are widely believed to be the key to 
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the formation of aspherical structure observed so commonly in post-AGB (pAGB) objects. 
The aspherical structure has been characterized via specific descriptors in a comprehensive 
morphological classification system, developed for PPNs and extended to young PNs (Sahai 
et al. 2007,2011). The descriptors primarily cover the shapes of the extended nebula and 
the central region. Two major classes of pAGB objects show equatorially-flattened central 
structures. The first class of pAGB objects consists of PPNs and PNs; a large fraction of 
these harbor overdense, dusty equatorial waists that are quite large; e.g., in several PPNs, 
the waists can be seen to display sharp outer (radial) edges, with radii typically > 1000 AU 
(Sahai et al. 2007). In the second class, namely dpAGB objects, there is strong evidence for 
medium-sized (~50 AU) disks (e.g., van Winckel et al. 2008) and little or no nebulosity. The 
relatively small CO line- width in Y Gem, implying the presence of a gravitationally bound 
structure, rather than a typical AGB dusty outflow, suggests that it will most likely evolve 
into a dpAGB object, rather than a PPN. 

Ultraviolet spectroscopic observations of Y Gem using HST's Cosmic Origins Spectro- 
graph are needed to probe the nature of the UV emission - i.e., if it is composed solely of 
emission lines arising in an accretion disk, or if there is a substantial continuum contribution 
(e.g., from the disk itself or the companion's photosphere). Any emission lines detected can 
provide constraints on the physical parameters of the accretion disk such as the (electron) 
density, emission measure, and thus the emitting volume (e.g., Reimers & Cassatella 1985). 
Multi-epoch observations (either at UV or radio wavelengths) can constrain the size of the 
emitting region depending on the time-scale of the variability. In a type 2b mechanism, one 
might expect gyroresonance or gyrosynchrotron emission at radio wavelengths from Y Gem 
- e.g., Skinner & Brown (1994) conclude that a likely model for the 2-3.6 cm emission in 
the pre- main- sequence star, TTau S, is non-thermal gyrosynchrotron emission arising in a 
scaled-up solar-like flare. Observations of high-J CO lines such as J=6-5, 7-6, 8-7, and 9- 
8 with the Herschel Space Observatory are necessary to determine the temperature of the 
molecular gas in the disk. 
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Table 1: GALEX UV Photometry of Y Gem fc Mira 
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Fig. 1. — (a) The variable UV emission from YGem. Both uncorrected values (circles) 
and corrected ones (for non-linearity) (triangles) of the FUV (blue symbols) & NUV (red 
symbols) fluxes vs. epoch are shown (the 0.15-mag error bars on the uncorrected values are 
the maximum zero point error that the GALEX images have had since the early mission 
data was released), (b) The observed SED: photometry from GALEX (UV), USNO-B.l 
and GSC2.3.2 (B,V,R), 2MASS (J,H,K), AKARI-MID (9, 18/xm), AKARI-FIS (90 /im), 
IRAS (12, 25, 60 /im), is shown as (black circles), and the IRAS/LRS spectrum as a black 
curve. A low-gravity T e ff = 2800 K model atmosphere spectrum from Lejeune et al. (1997), 
corresponding to Y Gem's M 8III spectral type (red curve) has been fitted to the SED. 
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Fig. 2. — CO J=2-l line emission from YGem, observed with the Arizona Radio Observa- 
tory's 10-m SMT mm-wave telescope. The temperature scale is measured as T\, where the 
source's radiation temperature is defined as T R = T\jr\^ and r] b (=0.74) is the main-beam 
efficiency. 



